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1.  SCIENTIFIC  WORK  DONE 


1 . 1  Logistics 

I  left  England  on  April  23,  1987  and  traveled  to  Vick¬ 
sburg,  Miss,  and  I  used  the  last  week  of  April  to  locate  ac¬ 
commodations  and  obtain  a  vehicle.  Although  the  contract 
did  not  start  officially  until  May  1,  I  did  visit  the 
Hydraulics  Lab.  at  Waterways  Experiment  Station  (WES)  prior 
to  that  date,  in  order  to  start  the  processes  of  gaining 
security  clearance  and  finding  desk  space  in  the  Hydraulic 
Analysis  Branch. 

1.2  Professional  Contacts 


When  I  began  work  on  May  1  my  first  task  was  to  make 
contact  with  the  other  scientists  and  engineers  working  on 
the  DEC  project,  and  on  other  projects  involving  the  study 
of  bank  stability.  This  was  essential  in  order  that  I 
quickly  learn  the  current  status  of  DEC  related  research, 
and  that  I  should  be  in  a  position  to  maximize  the  benefits 
of  work  already  completed,  by  applying  the  results  in  my  own 
study . 

At  WES  I  met  with  and  talked  at  length  to:  Mr  Bob 
Brown,  Mr  John  Ingram  and  Mr  Terry  Waller  in  the  Hydraulic 
Design  Branch.  Each  of  these  individuals  is  involved  in  DEC 
or  DEC  related  work.  Mr  Brown  is  a  member  of  the  DEC  Task 
Committee.  The  DEC  related  work  at  WES  is  concerned  with 
sediment  transport  prediction  and  stream  modeling  and  so 
bank  stability  considerations  are  of  immediate  relevance  and 
my  findings  will  be  of  interest  to  the  WES  staff. 

At  the  Vicksburg  District  I  have  worked  closely  with  Mr 
Phil  Combes  and  Mr  David  Biedenharn  in  the  Hydraulics  Sec¬ 
tion.  Both  are  directly  involved  in  the  DEC  Project,  and  Mr 
Biedenharn  isa  Corps  technical  advisor  to  the  DEC  Task  Com¬ 
mittee.  They  have  provided  vital  guidance  on  the  selection 
of  specific  locations  for  the  analysis  of  bank  stability, 
the  types  and  quality  of  historical  data  available,  and  the 
schedule  for  undertaking  the  work  in  this  project.  I  have 
also  worked  with  Mr  Chuck  Mendrop  in  the  Geotechnical  Sec¬ 
tion  at  the  District.  He  is  supplying  data  on  the  engineer¬ 
ing  properties  of  bank  materials  in  the  DEC  watersheds  from 
borehole  records  at  specific  locations. 

I  have  also  talked  with  other  individuals  who  are  work¬ 
ing  as  AE’s  on  the  DEC  Project.  Dr  Chester  Watson  of  Water 
Engineering  Technology,  Fort  Collins,  Co.  has  been  using 
bank  analyses  which  I  developed  in  1981,  to  undertake  a 
geomorphic  analysis  of  Hotophia  Creek.  He  has  found  bank 
stability  to  be  the  single  most  important  factor  in  control¬ 
ling  stream  channel  changes,  and  has  found  my  analysis  to  be 
successful  in  predicting  and  explaining  bank  response  to 
changing  hydraulic  and  sedimentary  conditions.  Dr  Watson 
has  made  available  to  me  all  of  the  results  of  his  studies 
to  date. 


Simons,  Li  and  Associates  have  recently  completed  their 
study  of  Hickahala/Senotobia  Creek.  Their  report  has  been 
made  available  to  me.  They  too  used  my  bank  stability  ap¬ 
proach,  but  only  qualitatively.  They  found  it  useful  in  the 
area  just  downstream  of  a  headcut  but  generally  give  the  im¬ 
pression  in  their  report  that  further  development  of  the 
analysis  is  needed  -  especially  to  include  additional  modes 
of  failure  such  as  rotational  slips.  This  will  be  incor¬ 
porated  into  my  current  research. 

Northwest  Hydraulic  Consultants,  Seattle,  Wa.  are  just 
beginning  studies  of  Black  Creek  and  of  Long  Creek.  I  have 
spoken  to  Mr  Charie  Neill,  the  PI,  and  told  him  of  my  work 
on  bank  stability  -  an  area  where  Northwest  have  compara¬ 
tively  little  experience.  We  plan  to  liaise  closely  in  the 
coming  months  to  incorporate  my  analyses  into  their  ap¬ 
proach  . 

In  the  first  months  of  this  project  I  have  developed 
good  working  relationships  with  the  individuals  working  on 
related  research  and  implementation  of  the  DEC  scheme.  This 
is  a  vital  prerequisite  for  the  success  of  my  project. 


1.3  Development  of  Bank  Stability  Theory 

Most  of  the  theoretical  analysis  was  already  completed 
prior  to  my  being  awarded  this  contract,  and  was  presented 
in  the  original  proposal  in  appendices  B  and  C.  Since 
starting  on  the  project,  I  have  further  developed  the 
theoretical  analysis  and  practical  applications,  specifi¬ 
cally  to  account  for  conditions  found  in  the  Mississippi 
Bluff  Line  streams.  The  theory  and  applications  are  set  out 
in  detail  in  papers  accepted  for  publication  in  the  Journal 
of  Hydraulic  Engineering  of  the  ASCE  (  Osman  &  Thorne,  1988; 
Thorne  &  Osman,  1988  ).  These  papers  have  been  taken  from 
draft  submission  to  final  form  during  the  first  two  months 
of  this  project  and  the  support  of  the  U.S.  Government  is 
duly  acknowledged. 

The  computations  involved  in  the  analysis  are  rela¬ 
tively  straight  forward,  but  still  present  numerous  oppor¬ 
tunities  for  arithmetic  errors.  Also,  correct  interpreta¬ 
tion  of  the  output  of  the  stability  model  demands  some  prior 
knowledge  of  stability  theory.  Therefore,  I  have  written  a 
program  for  the  Hewlett-Packard  HP-41C  calculator  which 
prompts  for  the  required  data,  performs  the  analysis, 
presents  and  interprets  the  results,  and  directs  the  user  to 
the  next  computational  step.  This  program  reouires  no  prior 
knowledge  of  the  HP-41C  and  can  be  used  on  site  as  the  cal¬ 
culator  is  battery  powered.  The  user  manual  gives  complete 
instructions  for  the  analysis  of  steep  streambanks  subject 
to  erosion  and  bed  degradation.  A  copy  of  the  user  manual 
is  attached  in  Annex  1. 


Testing  and  calibration  of  the  theoretical  analysis 
requires  field  data  from  unstable  banks.  These  data  have 
been  supplied  by  Mr  Biedenharn  for  sites  on  Long  Creek,  and 
its  tributaries  Caney  Creek  and  Goodwin  Creek.  Resurvey 
cross-sections,  long  profiles,  and  channel  planforms  for 
1978,1979,  and  1985  have  already  been  supplied  and  data  for 
1986  will  be  supplied  soon. 

2  RESEARCH  PLANS 


The  data  from  Long  Creek  will  be  used  to  test  predict¬ 
ions  of  bank  stability  and  channel  change  based  on  the 
theoretical  analysis.  Based  on  these  results  the  model  will 
be  modified  and  developed  as  necessary  to  produce  reliable 
predictions  of  bank  stability  on  Long  Creek.  The  results 
will  be  communicated  to  the  relevant  people  at  Northwest 
Hydraulic  Consultants,  for  incorporation  into  their  analysis 
of  the  Long  Creek  watershed. 

Next  it  is  planned  to  generalize  the  model  to  make  it 
applicable  to  other  watersheds  in  the  DEC  scheme  and 
throughout  Mississippi.  This  will  require  further  data  for 
testing  and  development. 

A  concurrent  aim  is  to  incorporate  the  bank  stability 
analysis  into  existing  Corps  models  for  aggradation  and 
degradation  in  streams.  This  major  undertaking  would 
require  considerable  effort  from  scientists  at  WES  in 
cooperation  with  myself,  but  could  greatly  improve  the  power 
and  utility  of  the  Corps  models.  I  have  put  this  proposal 
to  the  District  and  had  a  favorable  response  -  they  feel 
this  to  be  a  worthwhile  undertaking.  The  next  step  is  to 
talk  to  the  WES  personnel,  get  their  reaction  and  decide  how 
best  to  proceed.  This  is  an  immediate  priority. 


3  ADMINISTRATIVE  ACTIONS 


There  have  been  no  personnel  changes  or  important 
changes  of  an  administrative  nature  during  this  period.  The 
project  has  proceeded  as  set  out  in  the  original  proposal. 
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CALCULATOR  PROGRAM  FOR  ANALYSIS  OF  STREAMBANK  STABILITY 


Colin  R.  Thorne  Visiting  Scientist,  Vicksburg  District, 

US  Army  Corps  of  Engineers  (on  leave 
from  Queen  Mary  College,  Univ.  London, 
LONDON  El  4NS ,  England) 


1  INTRODUCTION 

Streambank  retreat  usually  occurs  by  a  combination  of  flow 
erosion  and  mass  failure.  Bank  failures  occur  when  erosion 
of  the  bank  and  the  bed  adjacent  to  the  bank  increase  the 
height  and  angle  of  the  bank  so  that  it  reaches  a  condition 
of  limiting  stability.  The  mechanism  of  failure  depends  on 
the  geometry  of  the  bank  and  engineering  properties  of  the 
bank  material  (Thorne , 1982 )  . 

Eroding  banks  are  usually  steep  and  commonly  fail  by  a 
slab-type  mechanism,  where  a  block  of  soil  topples  forward 
into  the  channel  (Fig.l).  The  weakening  effects  of  tension 
cracks  between  the  block  and  the  bank  can  be  important  in 
triggering  failure  and  should  be  accounted  for  when  analys¬ 
ing  the  stability  of  this  type  of  bank. 


Figure  1.  Slab  failure  of  a  streambank  in  Northwest  Missis¬ 
sippi  . 

Failures  usually  occur  during  "worst  case"  conditions, 
when  the  strength  of  the  bank  materials  is  minimised  and  the 
weight  is  maximised  due  to  high  moisture  levels.  Such  con¬ 
ditions  are  associated  with  periods  of  prolonged  rainfall, 
snowmelt,  and  drawdown  following  high  flow  stages  in  the 
channel.  Banks  which  are  stable  under  "average"  values  for 


1 


soil  properties,  but  which  are  unstable  for  "worst  case" 
conditions  are  at  risk  of  failure.  That  means  that  their 
stability  cannot  be  relied  upon  and  they  may  be  expected  to 
fail  sometime  in  the  near  future. 

The  analysis  of  the  stability  of  streambanks  with 
respect  to  slab  failure  was  undertaken  by  Thorne,  Murphey 
and  Little  (1981).  Recntly  Osman  and  Thorne  (1988)have 
developed  an  improved  approach  to  this  problem.  The  method 
uses  the  resolution  of  driving  and  resisting  forces  in 
static  equilibrium  on  the  most  critical  potential  failure 
plane  to  derive  a  factor  of  safety  ( FS )  for  a  bank  with 
respect  to  slab  failure.  A  value  of  FS  greater  than  one  in¬ 
dicates  stability,  and  equal  to  one  indicates  the  critical 
condition  with  the  bank  on  the  point  of  failure.  Values  less 
than  one  suggest  that  the  bank  is  unstable  and  should  have 
failed  already. 

This  manual  provides  a  calculator  program  for  the 
Hewlett-Packard  HP-41CV  or  CX  calculator,  which  may  be  used 
on-site  for  application  of  the  Osman-Thorne  analysis.  The 
HP-4 1C  may  be  used  provided  that  it  has  a  quad  memory  module 
f itted. 

This  Users  Guide  explains  how  the  program  is  used  to 
calculate  the  factor  of  safety  for  an  existing  bank,  and  how 
to  find  the  amounts  of  flow  erosion  and  bed  degradation 
which  would  cause  the  bank  to  fail.  The  bank  height,  ten¬ 
sion  crack  depth,  failure  plane  angle,  block  width,  and 
block  volume  for  slab  failure  are  computed.  These 
parameters  define  the  amount  of  top  bank  retreat,  the  volume 
of  slumped  material  input  to  the  channel,  and  the  new  bank 
geometry  after  failure. 

No  prior  experience  with  the  HP-41  is  necessary. 


2  DATA  REQUIREMENTS 

The  folowing  data  are  needed  to  perform  the  analysis: 


Variable 

Total  Bank  Height 
Upper  Bank  Height 
Bank  Slope  Angle 
Specific  Weight 
Effective  Friction  Angle 
Effective  Cohesion 
Tension  Crack  Index 


Symbol 

Units 

H 

ft 

HI 

ft 

I 

Degrees 

Sp . Wt . 

lb/ft 1 

phi 

Degrees 

c 

lb/ft 2 

K 

- 

The  tension  crack  index  is  defined  by: 

K  =  Y/H 
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where  Y  =  tension  crack  depth  (ft)  and  H  =  total  bank 
height.  Experience  shows  that  crack  depth  is  usually  limited 
to  less  than  half  the  bank  height.  If  there  are  no  data  on 
the  depth  of  cracking  at  the  site  being  investigated,  K  = 

0.5  may  be  used.  The  variables  are  shown  in  Figure  2. 


'JJ 
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Figure  2.  Definition  diagram  for  variables  in  slab  failure 
analysis:  a)  Initial  geometry;  b)  After  erosion  and  bed 
degradation  to  the  critical  case. 


3  INSTRUCTIONS 

First  check  if  the  program  is  already  in  the  cal¬ 
culator.  To  do  this  key  in: 

XEQ  ALPHA  USFS  ALPHA 

if  the  calculator  responds  with: 

NONEXISTENT 

then  you  mu3t  load  the  program  as  described  in  section  3.1. 
If  the  program  is  initiated  and  displays: 
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COLIN  THORNE-BANK  ANALYSIS 


then  you  may  go  straight  to  section  3.2. 


3.1  Loading  the  Program 

3.1.1  Manual  Loading 

1 .  Turn  on  the  calculator 

2.  Key  in  SHIFT  GTO  .  . 

3.  Put  the  calculator  in  Program  mode 

4.  Key  in  the  program  exactly  as  listed  in  appendix  4. 


Note , 

i)  All  statements  followed  by  either  ARCL  or  PROMPT  are 
entered  in  ALPHA  mode.  To  select  this  press  the  ALPHA  key. 
To  complete  entry  press  ALPHA  again. 

ii)  Commands  not  found  on  the  keyboard  (for  example 
PROMPT,  PSE  )  are  entered  by  keying; 

XEQ  ALPHA  command  ALPHA 

for  example  to  key  in  the  command  "PROMPT"  the  sequence  is; 
XEQ  ALPHA  PROMPT  ALPHA 

iii)  The  symbol,  -  ,  means  APPEND  and  is  keyed  in  by; 

ALPHA  SHIFT  K  - ALPHA 

for  example  to  key  in  the  units  that  follow  an  output  vari¬ 
able  the  sequence  is; 

ALPHA  SHIFT  K  FT  ALPHA 

iv)  The  SHIFT  key  is  yellow. 


3.1.2  Magnetic  card  loading 

If  you  have  the  program  on  magnetic  cards  and  have  a  card 
reader,  load  the  program  as  follows: 

1.  Turn  on  the  calculator 

2.  Key  in  GTO  .  . 

3.  Ensure  the  calculator  is  NOT  in  Program  mode 

4 .  Insert  the  cards  into  the  card  reader 

5.  When  all  cards  have  been  accepted,  key  in  GTO  .. 

The  program  is  now  ready  to  run. 
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3.2  Running  the  Program 

After  loading  the  program  (and  checking  for  keying  errors  if 
it  was  keyed  in  manually  )  the  program  is  initiated  by 
keying : 

XEQ  ALPHA  USFS  ALPHA 

After  that  it  is  only  necessary  to  press  the  RUN/STOP  (R/S) 
key  to  proceed,  entering  data  as  prompted  by  the  calculator 
and  noting  output  data  on  the  record  sheet.  Table  1.  Figure 
3  shows  the  program  flow  chart.  The  complete  user  instruc¬ 
tions  are  listed  in  section  3.2.1. 


3.2. 

1  User  Instructions 

STEP 

INSTRUCTION 

INPUT 

FUNCTION 

DISPLAY 

1 

Initialize  program 

XEQ 

USFS 

COLIN  THORNE- 
BANK  ANALYSIS 
STEP  NO.  1 . 

H  ?  FT 

2 

Input  total  height 

H 

R/S 

HI  ?  FT 

3 

Input  upper  height 

HI 

R/S 

I  ?  DEG 

4 

Input  bank  angle 

I 

R/S 

SP.WT  ?  LB/FT3 

5 

Input  specific  wt . 

Sp.Wt 

R/S 

PHI  ?  DEG 

6 

Input  friction  angle 

phi 

R/S 

c  ?  LB/FT2 

7 

Input  cohesion 

c 

R/S 

K  ? 

8 

Input  crack  index 

K 

R/S 

FS  = 

9 

Displays  Factor  of 
Safety  for  initial 
bank  geometry 

R/S 

A: 

BANK  STABLE 

B: 

H  =  CRITICAL 

10 

Displays  A  or  B 

If  A,  continue 

If  B,  go  to  step  19 

11 

Indicates  bank  is 
below  critical 
height, No  failure 
will  occur  without 
erosion  and/or 
degradation 

R/S 

STEP  NO. _ 

o 


12 

Displays  number  of 
current  step 

R/S 

DW  ?  FT 

13 

Input  lateral 

erosion  distance  DW 

R/S 

DZ  ?  FT 

14 

Input  degradation 
distance  DZ 

R/S 

H=  FT 

15 

Displays  new  H 

R/S 

HI  = _ FT 

16 

Displays  new  HI 

R/S 

FS  = 

17 

Displays  new  FS 

R/S 

A: 

BANK  STABLE 

B: 

H  =  CRITICAL 

C: 

BANK  TOO  HIGH 

18 

Displays  A,B,or  C 

If  A,  go  to  step  11 

If  B,  go  to  step  19 

If  C,  go  to  step  27 

19 

Indicates  bank  is  on 
point  of  failure 

R/S 

H  CRIT  = _ FT 

20 

Displays  critical  H 

R/S 

H1CRIT  = _ FT 

21 

Displays  critical  HI 

R/S 

Y  CRIT  = _ FT 

22 

Displays  critical  Y 

R/S 

I  CRIT= _ DEG 

23 

Displays  critical  I 

R/S 

BETA  = _ DEG 

24 

Displays  failure 
plane  angle,  BETA 

R/S 

BW  = _ FT 

25 

Displays  failure 
block  width,  BW 

R/S 

VB  = _ FT  3 /FT 

26 

Displays  failure 
block  volume  per 
per  unit  channel 
length,  VB 

R/S 

RUN  COMPLETE 

Program  returns  to 
step  1 

27 

Indicates  DW  &/or  DZ 
were  too  large  so 
bank  geometry  is 
unstable  and  could 
not  exist 

R/S 

REDUCE  DW  +  DZ 

6 


28 


Indicates  reduced 
DW  &/or  DZ  values 
must  be  used  and 

step  repeated  R/S  STEP  NO. 
Returns  to  step  12 


Table  1  Results  Table 


Stream  Name 

Site  Name 

Date 

Bank 

Initial  Data 

H  =  ft 

I  =  degrees 

phi  =  degrees 

HI  =  ft 

Sp.Wt  =  lb/ ft 3 

c  =  lb/ft* 

K  = 

Step 

DW 

(ft) 

DZ 

(ft) 

H 

(ft) 

H 

(ft) 

I 

(deg) 

FS 

COMMENTS 

i 

.  ■  _ 

Critical  Bank  Geometry  Parameters 

H  CRIT  =  ft 

I  CRIT  =  deg 

BW  =  ft 

H1CRIT  =  ft 

BETA  =  deg 

VB  =  ft’/ft 

BANK  ANALYSIS 


Figure  3.  Flow  Chart  for  Bank  Stability  Analysis 
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4  EXAMPLE 


4.1  Program  Steps 

KEYSTROKE  DISPLAY 


XEQ  ALPHA  USES  ALPHA 


15  R/S 
15  R/S 
70  R/S 
100  R/S 
12  R/S 
400  R/S 
0.5  R/S 
R/S 
R/S 
R/S 

0.6  R/S 

1.5  R/S 

R/S 

R/S 

R/S 

R/S 

R/S 

0.6  R/S 

1.2  R/S 

R/S 

R/S 

R/S 

R/S 

R/S 

R/S 

0.5  R/S 

1.1  R/S 

R/S 

R/S 

R/S 

R/S 

R/S 

R/S 

R/S 

R/S 

R/S 

R/S 

R/S 


COLIN  THORNE- 
BANK  ANALYSIS 
STEP  NO.  1. 

H?  FT 
HI?  FT 
I?  DEG 

SP.WT?  LB/FT3 
PHI?  DEG 
C?  LB/FT2 
K? 

FS  =  1.30 
BANK  STABLE 
STEP  NO .  2 . 

DW?  FT 
DZ?  FT 

H  =  16.50  FT 
Hl=  13.35  FT 
FS  =  1.11 
BANK  STABLE 
STEP  NO.  3. 

DW?  FT 
DZ?  FT 

H  =• 17.70  FT 
Hl=  11.70  FT 
FS  =  0.98 
BANK  TOO  HIGH 
REDUCE  DW  +  DZ 
STEP  NO.  3. 

DW?  FT 
DZ?  FT 

H  =  17.60  FT 
Hl=  11.98  FT 
FS  =  1.00 
H  =  CRITICAL 
H  CRIT  =  17.6FT 
H1CRIT  =  12. OFT 
Y  CRIT  =  8.8  FT 
I  CRIT  =  70  DEG 
BETA  =  44.7  DEG 
BW  =  4.54  FT 
VB  =  9 1 . 4FT3/FT 
RUN  COMPLETE 
COLIN  THORNE- 
BANK  STABILITY 
STEP  NO .  1  . 

H  ?  FT 
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4.2  Example  Results  Table 


Stream  Name  No  Name  Creek 

Site  Name  No  Place 

Date  6/18/87 

Bank  Left 

Initial  Data 


H  =  15.0  ft 

I  =  70  degrees 

Phi  =  12  degrees 

HI  =  15.0  ft 

Sp.Wt  =  100  lb/ f  t 3 

c  =  400  lb/ft3 

K  =0.5 

STEP 


DW 

(ft) 


DZ 

(ft) 


H 

(ft) 


HI 

(ft) 


I 

(deg) 


FS 


COMMENTS 


15.00 


15.00 


70 


1 . 30 


Stable 


0.6 


1 . 5 


16.50 


13.35 


70 


1.11 


Stable 


0.6 


1.2 


17.70 


11.70 


70 


0.98 


Too  High 


0.5 


1 . 1 


17.60 


11.98 


70 


1.00 


Critical 


Critical  Geometry  Parameters 


H  CRIT  =  17.6ft 


I  CRIT  =  70 


deg 


BW  =  4.54  ft 


H1CRIT  =  12.0ft 


BETA 


=  44.7  deg 


VB  =  91.4  ft’/ft 


Y  CRIT  =  8.8ft 
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APPENDIX  1 


SYMBOLS 


BETA 

c 

DW 

DZ 

FS 

H 

HI 

I 

K 

phi 

Sp.Wt 

VB 

Y 


Angle  between  the  failure  plane  and  the  horizon¬ 
tal.  Becomes  the  new  bank  angle  after  slab 
failure . 

Effective  cohesion  of  bank  material. 

Increment  of  lateral  bank  erosion  by  the  flow. 

Increment  of  bed  degradation  by  the  flow. 

Factor  of  safety  with  respect  to  slab  failure. 

Total  bank  height  above  the  bed. 

Upper  bank  height.  Equals  total  bank  height  if 
DW  and  DZ  are  zero.  Equals  total  bank  height  if 
bank  is  vertical. 

Angle  between  the  bank  surface  and  the  horizontal. 

Tension  crack  index.  K  =  Y/H. 

Effective  angle  of  internal  friction  of  bank 
material . 

Specific  weight  of  bank  material. 

Volume  of  failure  block  per  unit  channel  length. 
Depth  of  tension  cracking. 


DISPLAY 

BANK  UNSTABLE  : 
CHECK  DATA 


BANK  VERTICAL 
DUE  TO  EROSION 


HI > H  BAD  DATA 


W 1  <  0  ,  REDUCE  DW 


I >90,  BAD  DATA 


APPENDIX  2  -  MESSAGES  AND  ERRORS 
MEANING 


The  initial  bank  geometry  has  a 
factor  of  safety  less  than  1.  If 
"worst  case"  values  are  being  used 
then  this  means  that  the  bank  is 
at  risk  of  failure.  If  "average" 
values  are  being  used  then  the 
input  data  are  suspect  and  should 
be  checked. 

Keying  R/S  returns  the  program  to 
Step  1. 

Lateral  erosion  reduces  the  height 
of  the  upper  bank,  HI.  If  HI  goes 
to  zero  before  failure,  the  sloping 
portion  of  .the  bank  is  eliminated 
and  the  bank  is  made  vertical.  In 
this  case  the  program  makes  I  =  90 
and  HI  =  H.  Further  lateral  erosion 
produces  retreat  of  the  whole  bank 
and  does  not  affect  mass  stability. 
Therefore,  no  further  DW  values  are 
requested.  Vertical  banks  may  fail 
due  to  bed  degradation  represented 
by  further  DZ  increments. 

The  upper  bank  height  input  is 
greater  than  the  total  bank  height. 
This  is  impossible  and  the  program 
will  not  accept  H1>H.  Initially, 
if  DW  and  DZ  are  zero,  then  HI  =  H. 
Keying  R/S  returns  the  program  to 
Step  1. 

The  input  value  of  DW  would  reduce 
HI  below  zero  which  is  impossible. 
H<=H1>0  by  definition.  If  HI  goes 
to  zero  then  the  whole  bank  is 
vertical  and  H  =  HI.  This  is  taken 
into  account  in  the  program. 

Keying  R/S  returns  the  program  to 
Step  12  and  prompts  for  a  reduced 
DW  value. 

Input  value  of  the  bank  angle  I  is 
greater  than  90,  inferring  either 
data  error  or  an  cantilevered  bank. 
Cantilever  stability  can  simply 
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be  assessed  using  the  Thorne-Tovey 
(1981)  analysis.  The  slab  analysis 
used  here  would  be  inapplicable. 
Keying  R/S  returns  the  program  to 
Step  6  and  prompts  for  another 
value  of  I. 

K>1  BAD  DATA  Input  value  of  tension  crack  index 

K  is  greater  than  1,  which  is 
impossible.  Check  the  value  of  K. 
Keying  R/S  returns  the  program  to 
Step  8  and  prompts  for  another 
value  of  K. 

WARNING  :  I <60  Input  value  of  bank  angle  I  is  less 

60  degrees. Such  banks  often  fail  by 
rotational  slip,  not  slab  failure 
and  this  analysis  may  overestimate 
the  factor  of  safety.  Program 
continues,  but  results  should  be 
used  with  caution  for  I<60  degrees. 
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APPENDIX  4 


PROGRAM  LISTING 
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the  field  have  been  made  for  reconnaissance  and  site 
selection  for  application  and  detailed  testing  of  the 
method . 

In  July  I  traveled  to  Oxford,  Mississippi  to  visit  the 
L'SDA  Sedimentation  Laboratory  and  make  contact  with  the 
staff  there  who  are  involved  in  the  DEC  scheme.  Accompanied 
by  David  Siedenharn,  Charlie  Montague,  and  John  Smith  of  the 
LMVD ,  I  inspected  sites  of  severe  bank  erosion  on  Long, 
Marcum  and  Caney  Creeks,  and  looked  at  recently  constructed 
g rade-oont ro 1  structures  on  Long  and  Caney  Creeks.  Detailed 
data  collection  was  impossible  because  of  the  dense  growth 
of  vegetation  in  the  channel  and  on  the  bank  tops, 
especially  of  the  climbing  vine  "Kudzu" .  Although  this  vine 
lues  little  to  enhance  bank  stability,  it  does  limit  access 
to  the  banks  and  hides  areas  of  potential  instability  from 
v lew  . 

In  early  October,  I  traveled  to  the  Red  River  in 
Arkansas  and  Louisiana  with  David  Biedenharn,  Charlie  Little 
and  Freddie  Pinkard  of  the  LMVD.  The  purpose  of  this  trip 
was  to  examine  the  banks  of  the  Red  River  and  determine  the 
processes  and  mechanisms  of  failure.  It  is  possible  that 
the  Red  might  be  used  to  regionalize  the  results  of  the  bank 
stability  studies  in  the  DEC  watersheds. 

i.J  Professional  Contacts 

The  results  to  date  of  this  study  have  been 
communicated  to  interested  parties  at  WES  and  the  LMVD. 
Also,  liaison  with  the  AE’s  has  continued.  Specifically, 
copies  of  the  first  period  report  have  been  sent  to  Simons, 
Li  and  Associates  (SLA)  and  to  Water  Engineering  Technology 
(WET)  at  Fort  Collins.  The  project  has  also  been  discussed 
at  length  with  Dr  Mike  Harvey  of  WET.  No  feedback  has  been 
forth  coming  from  SLA.  Northwest  Hydraulics  have  withdrawn 
from  consideration  to  be  the  AE  for  Long  Creek,  and  no  .lE 
has  been  named  to  date. 

Data  on  channel  changes  and  geotechnical  properties  of 
the  bank  materials  has  been  supplied  by  the  LMVD  for  the 
years  1976,  l'J73,  1935  and  1986  for  various  reaches  of  Long 


Creek  and  its  tributaries,  Caney,  Goodwin,  Johnson  and 
Marcum  Creeks. 


1 . 3  Dose  Lopment  of  Ban k  Stability  Theory 

The  theoretical  and  applications  papers  describing  tr.e 
analysis  of  bank  stability  used  in  this  project  have  been 
scheduled  for  publication  in  the  February  1938  edition  of 
the  Journal  of  Hydraulic  Engineering  of  the  American  Soeietv 
of  Civil  Engineers  i Osman  and  Thorne,  1983;  Thorne  and 
Osman ,  1988a  )  . 

The  program  for  the  HP-llCV  calculator  has  been 
debugged  and  ex' ansively  tried  and  tested  for  accuracy. 
Minor  changes  have  been  made  to  the  program  to  elimmite 
hugs  and  to  improve  its  computational  efficiency.  V  new 
program  has  been  written,  to  analyse  subsequent  failures  due 
to  basal  clean-out  of  the  initial  failure,  and  further 
lateral  erosion  or  bed  degradation.  The  revised  and 
expanded  program  is  listed  in  Annex  1.  A  revised  us-rs’ 
manual  is  in  preparation. 

Cons  iderat  ion  has  been  given  to  the  role  of  batik 
stability  in  channel  evolution  towards  a  stable  aydrauiic 
geometry.  The  results  are  reported  in  a  paper  t,,  t-- 
presented  at  the  International  Conference  on  River  Regime  to 
bo  held  at  Wallingford,  UK  next  May  (Thorne  and  Osman, 
1933b).  Also,  some  thought  has  been  given  to  the  influence 
of  bank  vegetation  on  channel  geometry  and  particularly 
width.  The  results  are  reported  in  a  companion  paper  for 
the  same  conference  (Thorne,  et  al.  1988) 

1.1  Testing  the  Bank  Stability  Analysis 

On  the  basis  of  the  field  trip  to  Long  Creek  in  July,  a 
reach  around  structure  si,  at  range  300+00  upstream  of  the 
confluence  with  Johnson  Creek,  was  selectee  for  initial 
application  and  testing  of  the  bank  stability  analysis. 


Table  1 

.  Soil  parameters  for 

Long  Creek 

around  St 

rue turn  -i 

Soil 

percent  of 

Friction 

Cohes i on 

Ln  i  t 

Tensile 

Unit 

Bank  Height 

Ang  ie 

We l gh t 

Strength 

(%) 

(  ■>) 

( lb/ f t  * 

1  (lb/ft3) 

1 lb/ft3  ) 

PSA 

50 

12 

130 

123 

4  6 

MBI ( VP) 

20 

16 

317 

122 

1  06 

L  (OP) 

10 

1  1 

720 

132 

2  2  L 

Sands 

20 

20 

0 

1  1  1 

0 

46  • 


Bank  Averages 


16 


270 


130 


*  Only  upper  bank  tensile  strength  is  relevant  to 
calculation  of  tension  crack  depth  as  crack  depth  rarely 
exceeds  half  the  bank  height. 

The  data  on  soil  properties  supplied  by  the  L>1\  D  '-ere 
used  together  with  measurements  made  in  1330  on  the  banks  of 
Goodwin  and  Johnson  Creeks,  to  derive  representative  soil 
properties  for  the  four  soil  units  observed  in  the  banks  of 
Long  Creek  in  that  area.  The  mean  soil  properties  and 
weighted  averages  for  typical  banks  are  listed  in  Tabie  Z . 
Only  "worst  case"  values  are  listed  as  previous  experience 
shows  these  to  be  relevant  to  bank  stability  calculations 
iThorne  et  al.,  1981). 

The  data  presented  in  Table  2  were  used  to  estimate  the 
depth  of  tension  cracking  behind  the  banks  using  a  procedure 
developed  by  Thorne  et  al .  in  1981.  This  predicted  a 
maximum  crack  depth  of  5.16ft. 


Fig.l  Critical  Height  for  Slab-Type  Failure  of  Streambanks 
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The  soils  data  and  crack  depth  were  used  in  the  HP- 4  H  i 
program,  to  develop  curves  of  critical  bank  height,  critical 
degradation  depth  and  critical  lateral  erosion  distance  for 
the  study  reach.  Fig.  1  shows  the  critical  bank  height  as  a 
function  of  bank  angle.  Curves  based  on  the  earlier 
analyses  of  Chen  (i979)  and  of  Thorne  et  al .  (1981)  are 
plotted  for  comparison. 

Curves  to  predict  the  critical  degradation  depth  to 
instigate  mass  failures  under  gravity  for  banks  of  various 
initial  heights  and  angles  are  shown  in  rig.  2. 
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Critical  Degradation  Depth  for  Initial  Failure 


FrtiuuR E  Block 


Figs. 3  and  1  show  the  failure  block  width  and  block  volume 
per  unit  length  of  channel  respectively.  These  curves  could 
be  used  to  predict  the  channel's  width  response  tc. 
degradation  of  the  bed  and  estimate  the  volume  of  disturbed 
bank  debris  input  to  the  channel  as  a  result  of  bank 
instability. 


Fig. 3  Block  Width  for  Fig.  I  Block  Volume  for 
Initial  Failure  initial  Failure 

Following  the  initial  failure  due  to  degradation, 
further  failures  occur  when  the  flow  cleans-out  the  slump 
debris  and  erodes  the  bank.  Curves  of  critical  lateral 
erosion  distance  versus  initial  bank  angle  for  different 
initial  bank  heights  are  plotted  in  Fig.  5.  The  block  width 
is  in  this  case  equal  to  the  critical  erosion  distance.  The 
block  volume  is  obtained  from  rig  6.  These  curves  could  be 
used  to  predict  the  extent  of  widening  due  to  lateral 
erosion  and  basal  clean-out,  and  the  volume  of  sediment 
derived  from  bank  failures. 


way  to  determine  the  sequence  of  events  responsible  for 
bringing  the  bank  to  an  unstable  condition,  or  of  knowing 
the  geometry  at  the  time  of  failure.  Thus,  while  the  data 
are  of  use  for  general  testing  of  Fig.  1,  they  are  not 
sufficient  to  test  the  other  figures. 


bo  70  90  ,  SO  o 

/rgiTiA-u  AtJ&i-E  cnefigSE*) 

Fig. 6  Block  Volume  for  Subsequent  Failures  of  Long  Creek 
S  t  reambanks 

2  RESEARCH  PLANS 

To  overcome  the  limitations  of  the  data  supplied  bv 
LMVD ,  field  trips  will  be  made  to  Long  Creek  in  the  coming 
months,  to  make  detailed  observations  of  bank  geometry  just 
prior  to  and  .just  following  failure.  A  number  of  banks 
profiles  will  be  surveyed  in  areas  known  to  be  actively 


retreating,  and  the  staff  at  the  3ed imen tat i on  Lab'-rar  >ry 
will  keep  me  informed  of  high  flow  events  and  associated 
bank  erosion  leading  to  bank  failures,  so  that  banks  can  t.- 
re-surveyed  soon  after  failure. 

Regionalization  of  the  bank  analysis  using  data  from 
the  Red  River  is  distinctly  possible,  on  the  evidence  of 
this  month’s  visit  there.  This  possibility  will  be  pursued 
with  the  LMYD . 

Investigation  of  development  of  the  Corps  of  Engineers’ 
aggrada t lon/degrada t  ion  model,  HEC-6,  to  include  a  bank 
stability  criterion  will  continue. 


3  ADMINISTRATIVE  ACTIONS 

The  requirement  for  field  work  and  data  collection  has 
necessitated  that  a  field  assistant  be  available  to  me  to 
act  as  note  taker  and  surveying  rod-man.  Viso,  the  large 
volume  of  historical  channel  data  supplied  by  LM\ D  requires 
many  hours  of  work  on  data  reduction.  no  cost  to  the 
project,  my  graduate  student  at  Queen  Mary  College,  Lisa 
Cheadle,  will  help  me  with  these  tasks.  She  will  gain 
valuable  experience  and  exposure  to  channels  of  a  type 
unfamiliar  to  her,  which  will  benefit  her  in  her  stud:  -s 
toward  a  higher  degree  from  the  University  of  London. 
Therefore,  I  have  requested  of  the  European  Research  Offi.e 
that  Ms  Cheadle  be  added  to  the  project  as  a  temporary, 
unpaid,  field  assistant.  She  will  return  to  her  regular 
studies  in  January,  19S3. 
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48  STO  14 
♦9  2 

50  EHTERt 


51  1 

52  RCL  87 

53  - 

54  * 

55  RCL  06 

56  * 

57  RCL  04 
50  / 

59  RCL  02 

60  / 

61  STO  15 

62  RCL  14 

63  / 

64  XT2 

65  4 

66  ♦ 

67  S8RT 

68  RCL  15 

69  RCL  14 

70  > 

71  * 

72  2 

73  - 

74  STO  13 

75  RCL  29 

76  • 

77  STO  25 

78  FIX  2 

79  -r53=* 

30  ARCL  25 
3!  PRONP* 

32  1.895 

83  X<=Y" 

84  GTO  18 

85  RCL  25 

86  .995 

87  x>r 

38  GTO  19 

39  BEEP 

90  ”H1  =CR I T ICOL  - 

91  PROMPT 

92  FIX  2 

93  *HCRIT=" 

94  ARCL  01 

95  “EFT" 

96  PROMPT 

97  -mCRlT=- 

98  ARCL  02 

99  "EFT- 
100  PROMPT 


191  RCl  91 

182  RCl  87 

183  « 

104  STO  28 

105  "Y  CRIT=" 

106  ARCL  28 

107  "EFT- 

108  PROMPT 

109  "BETA=- 

110  ARCL  08 

111  "EBEG" 

112  PROMPT 

113  RCL  81 

114  RCL  82 

115  - 

116  RCL  08 

117  TAN 

118  / 

119  STO  38 

120  FIX  2 

121  "8H=" 

122  ARCL  36 

123  "EFT- 

124  PROMPT 

125  RCL  01 

126  Xf2 

127  RCL  02 

128  XT2 

129  - 
138  2 
13!  ’ 

132  RCL  83 

133  TAN 

134  • 

135  STO  31 

136  FIX  I 

137  -VB=* 

138  ARCL  3; 

139  -EFT3.--r- 

140  PRONPT 

14!  'RUN  CONFLET 

142  AVIEH 

143  PSE 

144  STOP 

145  GTO  -USFS- 
146HBL  18 

147  'BANK  STABLE 

148  PROMPT 
149*LBL  20 
150  1 
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151  ST*  24 

152  FIX  9 

201  RCL  01 

153  ‘STEP  NO.  ■ 

202  ♦ 

154  ORCL  24 

203  STO  01 

155  PROMPT 

204  FIX  2 

156  0 

205  “H=“ 

157  RCL  33 

206  ARCL  01 

158  -WP  FT- 

20?  -PFT- 

159  PROHPT 

208  PROMPT 

168  STO  2! 

209*LBL  24 

161  - 

218  *H1=- 

162  STO  33 

211  PRCL  82 

163  x>r> 

212  -FFT- 

164  GT0  2! 

213  PROMPT 

165  RCL  02 

214  CTO  1? 

166  RCL  32 

215M.8L  19 

167  - 

216  1 

168  RCL  21 

217  ST-  24 

169  RCL  08 

218  TONE  5 

170  T«K 

219  -BPNKTOOSTEE5 

171  * 

228  RVIEW 

172  X=0? 

221  PSE 

173  CTO  09 

222  -REDUCE  DN*53 

174  CLX 

223  PPOHP- 

175  CTO  21 

224*L8L  25 

176  TONE  5 

225  RCL  81 

177  *H1<0,REDUCEDH* 

226  RCL  22 

178  RRONPT 

179  RCL  21 

223  STO  ?! 

180  RCL  33 

229  RCl  >, 

181  + 

238  RCL  82 

182  STO  33 

231  * 

183  1 

232  STO  02 

184  ST-  24 

233  RCL  21 

185  CTO  20 

234  RCL  33 

186«L8L  21 

235  * 

187  RCL  21 

236  STO  33 

188  RCL  08 

23?  GTO  29 

189  TPM 

238  END 

190  * 

191  STO  16 

192  CHS 

193  RCL  02 

194  + 

195  STO  02 
196»L8l  22 

197  CLX 

198  •D2'>  FT 

199  PROMPT 

200  STO  22 
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DATE 

FILMED 


